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Abstract. Fractional velocity is defined as the limit of the difference quo¬ 
tient of the increments of a function and its argument raised to a fractional 
power. Fractional velocity can be suitable for characterizing singular behav¬ 
ior of derivatives of Holderian functions and non differentiable functions. The 
manuscript derives the product rules for fractional variation. Correspondence 
with integer-order derivatives is discussed. It is demonstrated that for Holder 
functions under certain conditions the product rules deviates from the Leibniz 
rule. This deviation is expressed by another quantity, fractional co-variation. 
Basic algebraic properties of the fractional co-variation are demonstrated. 
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1. Introduction 

Research on fractional derivatives and fractional calculus has long history m- 
Surprisingly, until the end of the 20*^ century only non-local definitions of fractional 
derivatives have been investigated. 

Applications to physical systems exhibiting fractal behavior have inspired the 
development of local definitions of fractional derivatives. Historically, the first non¬ 
integral order local definition has been proposed by Cherbit as the notion of a- 
velocity [9]. It was defined as a differential quotient of the increment of the particle 
displacement and the fractional power of the time interval between measurements. 
The concept implicitly implied Holder-continuity of the studied trajectory. The 
later definition of local fractional derivative due to Kolwankar and Gangal [12] 
was based on localization of Riemann-Liouville fractional derivatives, however this 
derivative is difficult to calculate for a wide range of functions. 

The correspondences between the integral and the quotient difference definitions 
have been recently investigated by Chen et al. [7]. On the other hand, Ben Adda 
and Cresson demonstrated equivalence between the a-velocity and the local frac¬ 
tional derivatives in the sense of Kolwankar and Gangal lais]. 

The present work uses the definition of a-fractional velocity of functions intro¬ 
duced by Cherebit. Results for a-velocity in this work are based on the previously 
introduced notion of fractional variation operators |18j . This manuscript demon¬ 
strates the product rule for a-velocity and establishes its relationship with the 
classical Leibniz’s product rule. 

The main results can be stated in the following pair of equations 
(x) = v^f {x) g{x) + v^g (x) f{x) ± [f,g]tix) , 
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where in turn v^f (x) designates the a-fractional velocity (or in short fractional 
velocity or a-velocity when discussing a particular order) of the function f{x) of 
order a and [/,(?]„ designates the fractional co-variation(s) of the functions f{x) 
and g{x) (of order a). Fractional velocity can be especially suitable for charac¬ 
terizing of singular behavior of derivatives of Hdlderian functions. Such functions 
are considered, for example, in quantum mechanical systems [I] or in the physical 
theory of scale relativity [mini Hz] where geodesic trajectories are considered to 
be non-differentiable. 

The present manuscript is organized as follows. Section[2]gives general definitions 
and notational conventions. Section introduces fractional variation operators 
and a-fractional velocities. Section |4] introduces fractional co-variation operators. 
Section [ 5 ] proves the main result and demonstrates the conditions under which the 
usual Leibniz rule holds. 

2. General definitions and notational conventions 

The term function denotes the mapping / : K 1 —>■ K (or in some cases C 1 —C). 
The notation f{x) is used to refer to the value of the function at the point x. The 
term operator denotes the mapping from one functional expression to another. The 
symbol C° denotes the class of functions, which are continuous. The symbol C” - 
the class of n-times differentiable functions where n G N. Square brackets are used 
for the arguments of operators, while round brackets are used for the arguments of 
functions. Dcym[f] denotes the domain of definition of the function fix). 

Definition 1. Let the parametrized difference operators acting on a function fix) 
be defined in the following way 

^tifKx) ■■= fix + e)-fix) , (1) 

KifKx) ■= fix) - fix- e ) , (2) 

^llfiix) ■= fix -be) - 2 fix) + fix - e) , (3) 

where e > 0. The first one we designate as forward difference operator, the sec¬ 
ond one as backward difference operator and the third one as order difference 
operator. 

Definition 2. Let IHI“ be the class of Holder 0^ functions of degree a, a G (0, 1), 
That is, we say that fix) is of class 11“ if V/(x) G 11“ there exist two positive 
constants C, <5 G M which for given x,y G Dom[f] such that for \x — y\ < 5 the 
following inequality holds 

\fix)-fiy)\<C\x-y\-, 

where \ ■ \ denotes the norm of the argument. Following Mallat and Hwang |14j for 
orders n > 1 the definition is extended in the following way: Holder class 
designates the class of 0 functions (of degree n + a) for which the inequality 

\fix) - fiy) - Puix -y)\<C\x- 2 /|”+“ , 
holds. Pni-) designates a real-valued polynomial of degree n G N 0/ the form 

n 

Pniz) = X] o.kZ^, where Poiz) = 0 and a G (0, 1]. 

k=l 

Remark 1. The polynomial Pnix) can he identified with the Taylor polynomial of 
order n of f{x) [T^ . 
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Under this definition we will focus mainly on functions for which 0 < a < 1. 
These functions will be further denoted by the term Holderian. 

Definition 3. Let the oscillation of a function f{z) in the interval [x, y\ be defined 
as 0 SC[ 2 ; y][/] := sup [/] — inf [/] for z € [x, y] or in alternative notation by 

Z ^ 

osc+[/] {x) := sup [f]{x) - inf [f]{x) 

e ^ 

in the interval [x^x -\- e]; or as 

osc7[/] {x) := sup [f]{x) - inf [f]{x) 

e « 

in the interval [x — e^x\. Let the oscillation about a point he given according to [21j ; 
osc[/](x) := lini sup [/] - inf [/] . 

[x-e,x+€] [x-£,x+e] 


3. Fractional variation and q!-fractional velocity of functions 


Fractional variation operators have been introduced in a previous work in the 
following way [18] : 


Definition 4. Let the Fractional Variation operators of order /3 he defined as 

[/] (,) := , (4, 

(5) 

where e > 0 and 0 < /3 < 1 are real parameters and f{x) is a function. The 
symbol vp~ [•] (•) will be called forward fractional variation operator of order f3 while 
[•] (•) will be called backward fractional variation operator. 


Remark 2. The following representation of the fractional variation operators is 
equivalent: 


vt [/] (^) 


A+[/](t) 


Ar[/](^) 


Remark 3. Following the notation of Cresson and Gref [ini[ii] we denote the 
usual forward and backward differential maps as: 

d; lf]M = IVlzIkzA . 

These maps, therefore, can be considered as particular cases of the Fractional vari¬ 
ation operators according to Definition\^ 


We define forward fractional velocity and backward fractional velocity, or in 
short fractional velocity or o-velocity when discussing a particular order, of order 
a in the following way: 
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Definition 5 (Fractional order velocity). Let the fractional velocity of fractional 
order a be defined as 

A+[/](x) 


!/(x) := lim 


e-)-0 


^ r [f]{x) 

vzf [X) := lim —— - 

e -)-0 e“ 


( 6 ) 

( 7 ) 


where e > 0 and 0 < a < 1 are real parameters and f(x) is function. 

Proposition 1 (Continuous variation). If f'{x) is continuous in the interval x £ 
[x, X -\- e] then 

v^f {x) = i lim e^~^f{x + e) , 
vtf (x) = ^ lim e^-^f{x - e) . 


Proof. By previous work 


□ 


The result is repeated for comparison with the expressions for fractional co¬ 
variation. 

Remark 4. In contrast to the restrictive treatment of the integral-ordered deriva¬ 
tives, we will not require equality of the forward and backward a-velocities. In such 
way it will be possible to handle symmetry breaking, which is important for physical 
applications, such as scale relativity where explicit symmetry breaking is admitted 

m- 

Indeed, if f{x) is defined in [x — e,x -f e] and the limit 

e-i-0 eP 

exists then v^f {x) ^ v^f (x). The proof of the statement is detailed in [18]. 

Example 1. We will calculate the forward and backward fractional variations of 
order 1/2 of f{x) = ^/x about x=0. 

^yx -I- e -I- y/x 
X -\- e — X e 


y/x -I- e -I- yfx y/x -\- C-\- y/x 
Then for a; = 0 [v^ i^) = 1- a similar manner 

a: m = = iff- ■ (v^+ 

X — X -\- e 


y/x - e -I- y/x 
e 


y/x — €-\- y/x y/x — €-\- yfx 
Then for a; = 0 u ^/2 [v^ ( 2 ^) = ~*- 

The same result can be computed using Proposition^^ In this case we have 
vV'^y/x (0) = 2 lim — r= = 1 

+ ^ ^ <=^0 2^0+7 

v^^'^y/x (0) = 2 lim — = —i 
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Theorem 1. Let f{x) he defined and finite in the interval [x, y\. Then 

\f{x) - f{y)\ < osc|^_y|[/] . 

Proof. The proof follows from the definition of the local oscillation operator. Let 
y = X + e. We will consider two cases. 

Case 1: Let f{x + e) > f{x). Since f{x + e) < sup [f]{x) and f{x) > inf [f]{x) 
then subtracting the inequalities leads to 

f{x + e)- f{x) <osc+[/](a;) . 

Case 2: Let f{x + e) < f{x). Since f{x) < sup [f]{x) and f{x + e) > inf [f]{x) 
then subtracting the inequalities leads to 

-fix + e) + fix) < osc+[/] ix) . 

Therefore, 

\fix + e)- fix)\ <osc+[/](a;) . 

□ 


Corollary 1. Let f{z) be monotone in [x,?/]. Then 

\fix) - fiy)\ = osC[^,y][/] . 

The proof follows immediately from the definitions and Theorem [T] 
Lemma 1 (Limit of the oscillatory quotient). 

osc+[/] (t) 


lim ■ 

e -)-0 

lim 

e-)-0 


osc, [/] (x) 


vff ix) 
V-fix) 


Proof. Let f{x) be monotone in [x,x + e] then the result follows trivially. Since e 
then is arbitrary there is an interval [x.,x + e'\, e' < e, such that fix) is monotone 
inside. The backward case can be proved in an identical manner. This completes 
the proof. □ 


Theorem 2 (Bounds of oscillation theorem). Let fix) G 111“, a € (0,1] in the 
interval [a:, a: + e]. Then 

0 < ce“ < osc+[/] (x) < C'e“ 
where c and C are positive real constants. 


Proof. There are two cases to be considered. 

Case 1: Let fix) be monotone. 

Then from the definition of Holder functions and Theorem [I] it follows 
that 


osc+[/] (x) < (76“ . 

On the other hand, osc+[/] (x) > B, where H is a non-negative constant. 
Dividing all inequalities by the positive amount e“ gives 


_B ^ osc+[/] (x) 
e“ “ e“ 


Since e can be made arbitrarily small then i? = c e“ in order to keep 
finite. 
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Case 2: Let f{x) be non-monotone. 

Then if we can split the interval in a collection of touching intervals 
where f{x) is monotone. 


Xn—X-\-e 


-I-e]\{a;-I-e} = |J [xi,x^+i) . 


Then, the same reasoning applies for every such interval and 
osc[:j._3,._^j][/] < CiAs" < Cie" . 

Therefore, 

osc^[/] (x) < max(C'i) e“ . 

Then the lower bound can be established by a similar reasoning as in the 
previous case. Then 

osC[a„_a„^p[/] > qAx" > 

where ki < AA-. Therefore, osc+[/] (x) > min(fci) e“. 

□ 

Proposition 2. Let f{x) S IHI“ in [x,x-|-e] then v^f {x) is unbounded if (3 > a, 
bounded if (3 = a and 0 if jd < a. 

Proof, [f] (x) < C e°‘~^. Therefore, the vanishing and the bounded cases 
follow. According to Theorem [2] 

osc+[/] (x) ^ 


Therefore, li f 3 > a the limit is unbounded. However by Lemma [T] 

vlf{x) 


= lim 


£->■0 e 

therefore, the fractional velocity is unbounded. 


□ 


The backward case is proven in a similar manner. 


4. Fractional co-variation of functions 

The product of two fractional variation operators will be designated as a new 
operator because of its utility in the expansion of products of a-velocities. 

Definition 6. Let the co-variation operators in the interval [x, x -|- e] be defined as 

:= [/] (x) la] (a^) , (8) 

where e > 0 and 0 < /? < 1 are real parameters and f(x) and g(x) are functions. 
When /3 = 1 the index will be omitted to keep notation simple. The symbol [•, will 
be called forward fractional co-variation operator while [•, will be called backward 
fractional co-variation operator. 

Remark 5. I can propose the following naming convention. The order ^ fractional 
co-variation to be named quadratic variation, the order ^ - ternary variation, the 
order -j - quaternary variation, etc. 
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Definition 7. We define the fractional co-variation of two functions f{x) and g{x) 
about a point x as the limit 

[/, g]^{x) := lim [/] {x) vfi^ [g] (x) (9) 

Proposition 3. The fractional co-variation is bi-linear. That is, without making 
discrimination about the direction of application for the numbers A, /i G M and the 
functions g{x),f{x), h(x) we have 

[A3 ,m/]/3 {x) = \g. [g,/]_a {x) , 

[a + f,h]p (x) = [g, hj^ (x) + [/, hj^ (x) . 

Proof. Both properties follows from the linearity of the fractional variation opera¬ 
tors. □ 

Proposition 4. The fractional co-variation is symmetric about its arguments. 
That is for functions g{x),f{x) without making discrimination about the direction 
of application we have 

[f.9\p{x) = [aJ\g{x) . 

The property follows from the commutativity of multiplication. 

Proposition 5. Let f{x) G and c is constant. Then without making discrimi¬ 
nation about the direction of application 

[/>c]^ (x) = 0 . 

The proposition is true since the difference operator yields zero from a constant. 



Figure 1. Fractional co-variation [f,gY^{x) for f{x) G and 
g{x) G ]HI“2 and /3 = 1. 


Theorem 3 (Bounds of co-variation). Let f{x) G and g{x) G Then 

lim[f,gY^(x) does not vanish only if ai-\-a 2 = fi. lim[/, is bounded only 

if ai -\- a 2 > fi- The result is illustrated in Fig. [IJ 
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Proof. Since the product of two Holder functions is a Holder function of the sum 
of the grades we have Af{x)Ag{x) G Therefore, by Theorem [2] it follows 

that 

C^C'e“i+“= < A+f{x)Atg{x) < C'iC 2 e“"+“" 
for certain constants C[,C 2 ,Ci,C 2 . Therefore, 

Taking the limit provides claimed result. The backward case follows from identical 
reasoning. □ 

Corollary 2. Let f(x) is uniformly C ^ m the interval [a;, a;+e] and g{x) £ ]HI“, a < 
1. Then without making diserimination about the direetion of application 

[/,5]/3 = 0 . 

Proof. By the last theorem for the exponents we have 1 + a > 1. This is satisfied 
for a > 0, therefore [/, g ]/3 vanishes. □ 


Computation of the fractional co-variation can be related to ordinary derivatives 
in the following way: 

Theorem 4 (Limiting case of Fractional co-variation). Let f{x),g{x) £ C^. Then 
lf’9]p = i^^^e^~^f'ix + e)Af[g]{x) + ^lime^-f^g'{x + e)At[f]{x) , 

^ p e—lO p e—>-0 

^ e)A7[g](a;) + ^ g'{x - e)Af[f]{x) . 

Proof. The proof follows from application of rHopital’s rule to the definition of 
fractional co-variation. For the forward case we have: 


[/,5]^ = lim 


€-)-0 
= lim 


{f{x + e)- f{x)) {g{x -k e) - g{x)) 

f{x + e)g[x -f e) - f{x)g{x -k e) - g{x)f{x -k e) -k g{x)f{x) 


In the last expression we treat e as an independent variable and a: as a parameter. 
Application of rHopital’s yields: 


[/> {f{x + e)g{x -k e) - f{x)g{x -k e) - g(x)f{x -k e) -k g{x)f{x))'^ 


A-t 


= 1™ + e) + + ^) 9 '{x -k e) - f{x)g'{x -k e) - g{x)f'{x + e)) 


e^O 


= ilimei ^f'{x + e)A+[g]{x) + ^\ime^ ^g'ix + e)At[f]{x) . 

P e-^O p e-^O 

In a similar manner, the backward co-variation can be computed as: 

If 1- f(^))i9ix-e)-gix)) 

[/, 9]p = l™,- -p - 

^ f{x - e)g{x - e) - f{x)g{x - e) - g{x)f{x - e) -k g{x)f{x) 

£-10 
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In the last expression we treat e as an independent variable and x as a parameter. 
Application of I’Hopital’s yields: 

[/> gTp = 1™ ifi^ - - e) - fix)g{x - e) - gix)f{x - e) + g(x)/(x))' 

= lim {-fix - e)g{x - e) - f{x - t)g'{x - e) + f{x)g'{x - e) + g{x)f{x - e)) 

£->0 p 

= 4 (a: - e)A-[g]{x) + 4 liin - e)A7[/](x) . 

P e-yO p £->-0 

□ 


There is also another reformulation possible. 

Corollary 3. Let /(x), g{x) G m the interval [x,x + e]. Then 

[f, 5 ]^ = ^ 1 ™ e /'(x + e) vf- [g] (x) + ^ ^ d'ix + e) [/] (x) , 

[/> 9]g = ^ 1™ e /'(x - e) Vp~ [g] (x) + ^ lim e ^'(x - e) [/] (x) . 

Proposition 6. For a function /(x) € in the interval [x,x + e] we have 

= I e /'(a; + e) [/] (x) , 

[/’ = I jiSo ■ 

Example 2 (Computation of fractional co-variation by Theorem|Tl). ITe will com¬ 
pute the a-velocity of the functions /(x) = ^/x and g(x) y/x for the degree 
a = 1/2. 




+ ^ 1 ^y/x + e-^ y/x y/x + e - ^/x 

1/2 — V ^ 77- 1 -V ^ ■ 




Vx + I 


3 3y^(x + e)2 


Eor X = 0 we have 


' e —V(l 




2 


33 


= 0 . 


For X ^ 0 since both /(x) and g{x) are continuous we have [f-,g]i /2 


= 0 . 


5. Product rule for fractional variation and o-velocity of 

HOLDERIAN FUNCTIONS 

In order to establish the main result we need to state two technical lemmas. The 
proofs of the lemmas are given only for completeness of presentation. 

Lemma 2 (P* Product Lemma). For a product of functions f{x)g{x) we have 

At [fg] (x) = A+ [/] (x) A+ [g] (x) + g(x)A+ [/] (x) + A+ [g] (x) f(x) . 
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Proof. Direct calculation shows that 
^^[fg] =f (a; + e) g(a; + e)-f(a;) g{x) 

= f (x + e) g (x + e) - f (x) g (x) + f (x + e) g(x) - f (x + e) g(x) 

= f (x + e) A+ [g] (x) + g(x)A+ [/] (x) 

= f (x + e) A+ [ 5 ] (x) + g(x)A+ [/] (x) + A+ [g] (x) f(x) - A+ [ 5 ] (x) f(x) 

= [/] ( 2 ;) [ 5 ] ( 2 ;) + g(x)A+f + A+ [g] (x) f(x) . 

□ 

Lemma 3 Product Lemma). 

A,” [fg] (x) = -A 7 [/] (x) A 7 [g] (x) + g(x)A 7 [/] (x) + A 7 [g] (x) f(x) . 
Proof. Direct calculation shows that 
^7 [f g] = f ( 2 ;) g (x) - f (x - e) g (x - e) 

= f (x) g (x) - f (x - e) g (x - e) + f (x - e) g(x) - f (x - e) g(x) 

= f (x - e) A 7 [g] (x) + g(x)A 7 [/] (x) 

= f (x - e) Aj [g] (x) + g(x)A 7 [/] (x) + A 7 [g] (x) f(x) - A 7 [g] (x) f(x) 
= [/] ( 2 ;) A 7 [g] (x) + g(x)A 7 [/] (x) + A 7 [g] (x) f(x) . 

□ 


Having established these results, the product rule for the fractional variation 
operators and a-velocity can be stated in the following Theorem: 

Theorem 5 (Product rule for Fractional variation). Let f{x) G 11“^ and g{x) G 
Then 


[/ 9 ] {x) = vf- [/] (x) g(x) + vf- [g] (x) /(x) + [/, g]^+ (x) , 
3 ] {x) = Vp~ [f] (x) g{x) + Vp~ [g] (x) /(x) - [/, g]^" (x) . 

Proof. We prove the forward case first. 

Forward case: By Lemma [5] it follows that 

A+[fg]( 2 :) = A+[f](x) A+[g](x) +g(x)A+[f](x) + A+[g](x)f(x) 
Dividing the last quantity to yields 


( 10 ) 

( 11 ) 


A+[fg](x) _ A+f(x)A+g(x) 


V [-C] (^) 9(x) + V ^ 3 ] (x) fix) 


gp gp 

= [/: ( 2 :) + [/] (x) g(x) + [g] (x) fix) , 

which completes the proof. 

Backward case: The backwards case follows from analogous reasoning. By 
Lemma [3] it follows that 

Ar[f g]( 2 ;) = -A 7 f(x) A 7 g(x) + g(x)Ai,f(x) + A 7 g(x) f(x) . 

Dividing the last quantity to yields 

A£[^ = [/] (^) 3(:g) + [ 5 ] (^) /(^) 

= -[f,9]p~ix) + V 0 ~ [/] ( 2 :) g(x) + [g] (x) /(x) , 
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( 12 ) 

(13) 


which completes the proof. □ 

Corollary 4 (Limiting case). For the functions f{x) and g{x) 

v+[f a] {x) = v^f {x)g{x) + vlg{x) f{x) + [f,g\p{x) , 

'^-[f a] {x) =vtf{x)g{x) + vtg{x)f{x)- [f,g\~p{x) . 

This corollary confirms a recent result by Tarasov [20] who affirms violation of 
the Leibniz’s rule for fractional derivatives of non-differentiable functions. Indeed, 
by Theorem [3] the quantity [/, is not identically zero. 

Proposition 7 (Square case). For the function f{x) we have 

v^f (x) = 2 vlf {x) f{x) + [/, f]^ (x) , 

vtf (x) = 2 vtf (x) f(x) - [/, f]^ (x) . 

Corollary 5 (Quotient case). For the functions f{x) and g{x) 

I a] {x) = ^ , 

Pfr, w N v^f(x)gix)-v^g{x)f(x) + {f,g]~ 

V- [fla] {.x) = -^ , 

where we assume that g(x) ^ 0. 

Proof. We will prove the forward case first. Since g(x) ^ 0 we have 
V%f (x) ^ + v% [1/g] (x) /(x) + [/, l/g\^ (x) . 

The fractional variation of the quotient function can be calculated as 


(14) 

(15) 

(16) 
(17) 


1 


1 


C, \y!i]{^) = T>\7PT7)-W) 

Taking the limit yields 

[a] (x) 


,e+ 


la] (x) 


< [ 1 / 5 ] (a;) = - lim 


a{x)g{x 

v+a(x) 


£^‘6 g{x)g{x + e) g'^{x) 
since by assumption g(x) ^ 0. Therefore, for the fractional co-variation it follows 


that 


and 


a ^^/2 [/] {x) [g] (x) 


[/, Va]t{x) = 

P g{x)g{xFe) g{x)g[x + e) 




£-1- 


respectively. 

For the backward case we have to not that 


Taking the limit yields 


[a]ix) 


g{x)a{x - e) ■ 
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since by assumption g{x) ^ 0. Therefore, for the fractional co-variation it follows 
that 

. , , _ ^'p/2 [/] (^) la] i^) _ [/,g]r 

^ g{x)g{x-e) g{x)g{x - e) 

and taking the limit finally yields 

respectively. □ 

The proof established also the following additional results: 

Proposition 8. For the functions f{x) and g{x) ^ 0 we have 

[fA/9\p{.x) = --^4^ • 

P g-^[x) 

Proposition 9. For the function g{x) ^ 0 we have 

/5 n / 1/ , ^±9{^) 

"±[V9lW = — 

Example 3 (Direct calculation of quadratic variation). We will calculate the qua¬ 
dratic variation of f{x) = ^/x. Then 

[yJxF e - i/i) [s/xF e + i/i) _ 


A+ f[x) = y/xT+l. — ^/x = 


X + € — X 


\/x F e -I- -v/x 
e 


\/x -I- e -I- sjx \/x -I- e -I- sjx 


Therefore 


Wx, = 


(xFtt-)' 


Therefore, 


^ [y/xF e -I- ^/x) 

Example 4 (Calculation of quadratic variation from Theorem |4|). 


^ 

[v/x, = li™, 2 ;; — ! {^x F e — \/x) = 1 — lim ■ 


e—>0 2 X F > 


£->■0 ^x F I 


Therefore, 


[y/x, v^]+ = I 


1, X = 0 , 

0, X > 0 . 


Combining the evaluation of the product rules leads to the following results 
Proposition 10. For functions f{x) and g{x) 
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Proposition 11. For functions f{x) and g{x) 

vt [/ g] {x) lim {g'{x - e)f{x) + g{x)f'{x - e)) - 

p e->0 

i Im {g\x - e)A7 [/] (x) + f{x - e)A7 [g] (x)) 

5.1. Relation to the Leibniz’s rule. Under the particular notation in Remark 
13] the product rule restated in the following way: 

Proposition 12. Under the notation used in this work, the product rule for ordi¬ 
nary differential quotients can be expressed by 

[/ g] {x) = d+ [/] (x) g{x) + d+ [g] (x) /(x) + [/, g]+ (x) , (18) 

dr [/ g] (x) = dr [/] (x) g{x) + dr [g] (x) /(x) - [/, 3]r (a:) . (19) 

Proof. The result has been stated in [TOl [11]. However, we can consider it as a 
specialization of the a-fractional velocity product rule (Theorem [5]) for the case 
a = 1. □ 

If we take the limit we arrive at the usual Leibniz rule for derivatives: 

Proposition 13. 

d"^ [fg] (t) = d+ [/] (x)g(x) +d+ [g] (x)/(x) + [f,g]^{x) (20) 

d" [fg] (x) = d“ [/] (x)g(x) +d" [g] (x)/(x) - [/,g]“(x) (21) 

For functions in with uniformly continuous derivatives taking the limit to zero 
yields the usual Leibniz rule since [f,g]~{x) = [/,g] + (x) = 0. 

6. Discussion 

The theory of local fractional derivatives is still immature and there are few 
established results with certainty n IS1I31 mu mis]. Historically, the first non¬ 
integral order local definition has been proposed by Cherbit [5] as the notion of 
” a-velocity” in the study of fractals as 

X(«) := lini + 

s^O S°‘ 

where a G (0,1] for the trajectory X(t) as a function of time. The concept 
implicitly implied continuity. Chen et al. used the term ’’difference-quotient” 
based local fractional derivative [7]. It should be clear that the terms a-velocity 
(Cherebit), ’’difference-quotient” based local fractional derivative (Chen et al.) and 
” a-derivative” (Ben Adda and Cresson) are equivalent. Regrettably, in the current 
literature there is still no unification of notations and terms, which provides space 
for confusion. Properties of this entity are qualitatively different from the proper¬ 
ties of ordinary derivatives and the integral fractional derivatives. Therefore, its 
study has its own merits. For example, it has very few non-zero or non-divergent 
values therefore, it has more of a supplementary character compared to integral 
derivatives. 
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